ITP-Budapest 560 



Grand Unification signal from ultrahigh energy cosmic rays? 
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The spectrum of ultrahigh energy (above « 10^ GeV) cosmic rays is consistent with the decay of 
GUT scale particles. The predicted mass is mx ~ IQ^ GeV, where h = 14.6li;?. 



PACS numbers: 98.70.Sa, 96.40.Pq, 13.87.Fh 



O 
O 



(N • 

> : 
^ . 
o ■ 

00 

o 

O 

o : 

i!^ • 
Oh; 

I 

Oh 



X 



The interaction of protons with photons of the cos- 
mic microwave background radiation (CMBR) predicts 
a sharp drop in the cosmic ray flux above the Greisen- 
Zatsepin-Kuzmin (GZK) cutoff around 5 • lO^^ eV 0. 
The available data shows no such drop. About 20 events 
above 10^° eV were observed by experiments such as 
Akeno Giant Air Shower Array (AGASA) Fly's Eye 
H, Haverah Park [g, Yakutsk |] and HiRes g. Future 
experiments, particularly Pierre Auger Q, will have a 
much higher statistics. 

Usually it is assumed that at these energies the galactic 
and extragalactic (EG) magnetic fields do not affect the 
orbit of the cosmic rays, thus they should point back to 
their origin within a few degrees. Though there are clus- 
tered events |^,^ the overall distribution is practically 
isotropic ||l0|| , which usually ought to be interpreted as a 
signature for EG origin. 

Since above the GZK energy the attenuation length of 
particles is a few tens of megaparsecs pl]-[T^ if an ultra- 
high energy cosmic ray (UHECR) is observed on Earth 
it must be produced in our vicinity (except for UHECR 
scenarios based on weakly interacting particles j_e.g. neu- 
trinos (l^). Sources of EG origin (e.g. AGN JM, topo- 
logical defects |17 or the local supercluster [|18[) should 
result in a GZK cutoff, which is in disagreement with ex- 
periments. It is generally believed that there is no 
conventional astrophysical explanation for the observed 
UHECR spectrum. 

An interesting idea suggested by refs. |po| , pT| is that 
superheavy particles (SP) as dark matter could be the 
source of UHECRs. (Note, that metastable relic SPs 
were proposed before the observation of UHECRs 
beyond the GZK cutoff.) In EG SPs were stud- 
ied. Ref. made a crucial observation and analyzed 
the decay of SPs concentrated in the halo of our galaxy. 
They used the modified leading logarithmic approxima- 
tion (MLLA) [p3| for ordinary QCD and for supersym- 
metric QCD [p4||. A good agreement of the EG spectrum 
with observations was noticed in p5[ |. Supersymmctric 
QCD is treated as the strong regime of the minimal su- 
persymmetric standard model (MSSM). To describe the 
decay spectrum more accurately HERWIG Monte-Carlo 
was used in QCD |^6|] and discussed in supersymmctric 
QCD HJH, resulting in mx ~ lO^^ Qgy and « lO^^ 
GeV for the SP mass in SM and in MSSM, respectively. 



SPs are very efficiently produced by the various mech- 
anisms at post inflatory epochs ||2^ . Note, that our anal- 
ysis of SP decay covers a much broader class of possible 
sources. Several non-conventional UHECR sources (e.g . 
EG long ordinary strings or galactic vortons 
monopole-antimonopole pairs connected by strings [p2[) 
produce the same UHECR spectra as decaying SPs. 

In this letter we study the scenario that the UHECRs 
are coming from decaying SPs and we determine the mass 
of this X particle mx by a detailed analysis of the ob- 
served UHECR spectrum. We discuss both possibilities 
that the UHECR protons are produced in the halo of our 
galaxy and that they are of EG origin and their propa- 
gation is affected by CMBR. Here we do not investigate 
how can they be of halo or EG origin, we just analyze 
their effect on the observed spectrum instead. We assume 
that the SP decays into two quarks (other decay modes 
would increase mx in our conclusion). After hadroniza- 
tion these quarks yield protons. The result is character- 
ized by the fragmentation function (FF) D{x, Q^) which 
gives the number of produced protons with momentum 
fraction x at energy scale Q. For the proton's FF at 
present accelerator energies we use ref. [Q. We evolve 
the FFs in ordinary [Q and in supersymmctric QCD 
to the energies of the SPs. This result can be combined 
with the prediction of the MLLA technique , which gives 
the initial spectrum of UHECRs at the energy mx ■ Al- 
together we study four different models: halo-SM, halo- 
MSSM, EG-SM and EG-MSSM. 

Ref. |§ showed that both AGASA and Fly's Eye data 
demonstrated a change of composition, a shift from heavy 
-iron- at lO^"^ eV to light -proton- at 10^^ eV. Thus the 
UHECRs are most likely to be dominated by protons and 
in our analysis we use them exclusively. 

The FF of the proton can be determined from present 
experimen ts pd^ . (Note, that QCD event generators e.g. 
HERWIG |37| predict the overall proton multiplicity cor- 
rectly, however they describe the large x region of the FF 
inaccurately.) The FFs at Qo energy scale are Di{x,Q'^), 
where i represents the different partons (quark/squark or 
gluon/gluino). The FFs can not be determined in pertur- 
bative QCD; however, their evolution in is governed 
by the DGLAP equations [|3: 
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u = 2d 


1.41 


0.402 


-0.860 


2.80 


s 


1.41 


4.08 


-0.0974 


4.99 


c 


2.9 


0.111 


-1.54 


2.21 


b 


9.46 


40.1 


0.742 


12.4 


t 


350 


1.11 


-2.05 


11.4 


g 


1.41 


0.740 


-0.770 


7.69 




1000 


0.82 


-2.15 


10.8 



TABLE I. The fragmentation functions of the different 
partons using the parametrization D{x) = Nx°'(l — x)^ at 
different energy scales (second column). 
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One can interpret Pji{z), the splitting function, as the 
probability density that a parton i produces a parton j 
with momentum fraction z. Analogous evolution equa- 
tions can be obtained by using coherent branching (angu- 
lar ordering) for the emitted gluons [ ^8|j2^ . We use this 
technique too. Results of direct Monte-Carlo jet simu- 
lations are also available (cf. [^). We include the un- 
certainties coming from the different choices of evolution 
in our final error estimates. We solve the DGLAP equa- 
tions numerically with the conventional QCD (SM case) 
splitting functions and with the supersymmetric (MSSM 
case) ones |Q. For the top and the MSSM partons we 
used the FFs of ref. Hi. While solving the DGLAP equa- 
tions each parton is included at its own threshold energy. 
Table | shows all the FFs we used. 

At small values of x, multiple soft gluon emission 
can be described by the MLLA |^^. This gives the 
shape of the total hadronic FF for soft particles (not 
distinguishing individual hadronic species) xF{x,Q^) oc 
exp [— ln(a;/a::m)2/(2cr2)] , which is peaked at Xm — 

^/A/Q with 2ct2 = Aln^/^(Q/A). According to || the 
values of A are ^/fJS/G and 1/6 for SM and MSSM, re- 
spectively. The MLLA describes the observed hadropro- 
duction quite accurately in the small x region |4C[ |. For 
large values of x the MLLA should not be used. We 
smoothly connect the solution for the FF obtained by 
the DGLAP equations and the MLLA result at a given 
Xc value. Our final result on mx is rather insensitive to 
the choice of Xc, the uncertainty is included in our error 
estimate. We also determined the FF of the pion. Fig. |l| 
shows the FF for the proton and pion at Q = 10^^ GeV 
in SM and MSSM. 

UHECR protons produced in the halo of our galaxy 
can propagate practically unaffected and the production 
spectrum should be compared with the observations. 




FIG. 1. The FFs averaged over the quark flavors at 
Q = 10^^ GeV for proton/pion in SM (solid/dotted hne) and 
in MSSM (dashed/dashed-dotted line) in the relevant x re- 
gion. To show both the small and large x behavior we change 
from logarithmic scale to linear at x = 0.01. 



Particles of EG origin and energies above w 5 • 10^^ eV 
loose a large fraction of their energies due to interactions 
with CMBR jlj. This effect can be quantitatively de- 
scribed by the function P{r, E, Ec), the probability that 
a proton created at a distance r with energy E arrives at 
Earth above the threshold energy Ec . This function 
has been calculated for a wide range of parameters in 
p2[ , which we use in the present calculation. The origi- 
nal UHECR spectrum is changed at least by two different 
ways: (a) there should be a steepening due to the GZK 
effect; (b) particles loosing their energy are accumulated 
just before the cutoff and produce a bump. We study the 
observed spectrum by assuming a uniform source distri- 
bution for UHECRs. 

Our analysis includes the published and the unpub- 
lished (from the www pages of the experiments on 
17/08/00) UHECR data of ||-|,§. Due to normahza- 
tion difficulties we did not use the Yakutsk |^ results. 
We also performed the analysis using the AGASA data 
only and found the same value (well within the error 
bars) for mx- Since the decay of SPs results in a non- 
neghgible flux for lower energies log(i5mm/eV) = 18.5 
is used as a lower end for the UHECR spectrum. Our 
results are insensitive to the definition of the upper end 
(the flux is extremely small there) for which we choose 
log(i?maa;/eV) = 26. As it is usual we divided each loga- 
rithmic unit into ten bins. The integrated flux gives the 
total number of events in a bin. The uncertainties of the 
measured energies are about 30% which is one bin. Us- 
ing a Monte-Carlo method we included this uncertainty 
in the final error estimates. The predicted number of 
events in a bin is given by 
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FIG. 2. The available UHECR data with their error bars 
and the best fit from a decaying SP. Note that there are no 
events above 3 x 10^" eV (shown by an arrow). Nevertheless 
the experiments are sensitive even in this region. Zero event 
does not mean zero flux, but a well defined upper bound for 
the flux (given by the Poisson distribution). Therefore the 
experimental value of the integrated flux is in the "hatched" 
region with 68% confidence level, ("hatching" is a set of in- 
dividual error bars; though most of them are too large to be 
depicted in full) Clearly, the error bars are large enough to be 
consistent with the SP decay. 
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where Ei is the lower bound of the i*'' energy bin. The 
first term describes the data below 10^^ eV according to 
1^ , where the SP decay gives negligible contribution. The 
second one corresponds to the spectrum of the decaying 
SPs. A and B are normalization factors. 

The expectation value for the number of events in a 
bin is given by eqn. (||) and it is Poisson distributed. 
To determine the most probable mx value we used 
the maximum-likelihood method by minimalizing the 
X^iA, B,mx) for Poisson distributed data 

26.0 

X'= 5] 2 [N{t)-No{t) + No{t) In {No{t)/Nm, (3) 

i=18.5 

where No{i) is the total number of observed events in the 
•t/i ^jjj Q^j. g^^jjjg procedure we have three parame- 
ters: A, J5 and mx- The minimum of the ~x^{A,B,mx) 
function is Xmm mxmin which is the most probable 
value for the mass, whereas x^{A' , B' , mx) = Xo(™^) — 
Xmin + 1 gives the one-sigma (68%) confidence interval 
for mx- Here A' , B' are defined in such a way that the 
X^iA, B, mx) function is minimalized in A and B at fixed 
mx- Fig. 1^ shows the measured UHECR spectrum and 
the best fit in the EG-MSSM scenario. The first bump 
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FIG. 3. The most probable values for the mass of the 
decaying ultra heavy dark matter with their error bars and 
the total values. Note that 21 bins contain nonzero number 
of events and eqn.(S) has 3 free parameters. 



of the fit represents particles produced at high energies 
and accumulated just above the GZK cutoff due to their 
energy losses. The bump at higher energy is a remnant 
of mx- In the halo models there is no GZK bump, so 
the relatively large x part of the FF moves to the bump 
around 5 x 10^^ GeV resulting in a much smaller mx than 
in the EG case. An interesting feature of the GZK effect 
is that the shape of the produced GZK bump is rather 
insensitive to the injected spectrum so the dependence of 
on the choice of the FF is small. The experimental 
data is far more accurately described by the GZK effect 
(dominant feature of the EG fit) than by the FF itself 
(dominant for halo scenarios). 

To determine the most probable value for the mass of 
the SP we studied 4 scenarios. Fig. ^ contains the Xmin 
values and the most probable masses with their errors for 
these scenarios. (The uncertainties coming from the FFs 
are included in our error estimates on mx-) 

The UHECR data favors the EG-MSSM scenario. The 
goodnesses of the fits for the halo models are far worse. 
The SM and MSSM cases do not differ significantly The 
most important message is that the masses of the best 
fits (EG cases) are compatible within the error bars with 
the MSSM gauge couphng unification GUT scale Q. 

The SP decay will also produce a huge number of pi- 
ous which will decay into photons. Our spectrum con- 
tains 94% of pions and 6% of protons. This Ti/p ratio is 
in agreement with which showed that for differ- 

ent classes of models mx ^ 10^^ GeV , which is the up- 
per boundary of our confidence intervals, the generated 
gamma spectrum is still consistent with the observational 
constraints. We performed the whole analysis including 
the pion produced 7-s in eqn. (^. The results agree 
with our results of Fig. ^ within errorbars, which is easy 
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to understand. For the EG case high energy 7-s domi- 
nate at energies where the observed flux is zero . For 
the halo case the agreement is resulted by the similarity 
(except normalization) between Dp and (cf. Fig. l|). 

In the near future the UHECR statistics will probably 
be increased by an order of magnitude 0|. Performing 
our analysis for such a statistics the uncertainty of mx 
was found to be reduced by two orders of magnitude. 

Since the decay time should be at least the age of the 
universe it might happen that such SPs overdose the 
universe. Due to the large mass of the SPs a single decay 
results in a large number of UHECRs, thus a relatively 
small number of SPs can describe the observations. We 
checked that in all of the four scenarios the minimum 
density required for the best-fit spectrum is more than 
ten orders of magnitude smaller than the critical one. 

Details will be presented in a subsequent paper . 

We thank B.A. Kniehl for providing us with the pro- 
ton's FF prior to its publication and F. Csikor for 
useful comments. This work was partially supported 
by Hungarian Science Foundation grants No. OTKA- 
T29803 /T22929-FKP-0128/1997. 
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